Optimizing the catching of atoms or molecules in two-dimensional traps 
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Single-photon cooling is a recently introduced method to cool atoms and molecules for which stan- 
dard methods might not be applicable. We numerically examine this method in a two-dimensional 
wedge trap as well as in a two-dimensional harmonic trap. An element of the method is a small op- 
tical box with "diodic" walls which moves slowly through the external potential and catches atoms 
irreversibly. We show that the cooling efficiency of the method can be improved by optimizing the 
trajectory of this box. 
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I. INTRODUCTION 

Recently, a new method has been introduced, called 
single-photon cooling 0, Q , which allows to cool atoms 
and molecules which cannot be handled in a standard 
way. The method is based on an atom diode or one- 
way barrier 0, |4j. An atom diode is a device which 
can be passed by the atom only in one direction whereas 
the atom is reflected if coming from the opposite direc- 
tion. Such a device has been studied theoretically [l-ll 
and also experimentally implemented as a realization of 
a Maxwell demon [9|, [lO( . 

The idea of cooling based on an atom diode is ex- 
plained in [1| and [7|. For illustration, we assume a sim- 
plified, idealized one-dimensional setting and an ensem- 
ble of non-interacting atoms which are moving classically 
in a harmonic potential trap, see Fig. [I] During the whole 
cooling process the atom diode is moved slowly through 
the trap from right to left with a constant velocity. In 
addition, we assume that the width of the diode is neg- 
ligible. An atom can cross the atom diode only from left 
to right and after such a crossing the atom is trapped on 
the right-hand side of the diode. 

For simplification let us look in the following at the pro- 
cess for a single atom of the ensemble, see Fig. [1] and also 
0. As the diode is moving much slower than the atom 
the diode captures every atom near its classical turning 
point, i.e. the atom has very low kinetic energy when 
it is caught. If the atom crosses the diode from left to 
right then this is an irreversible process, the atom cannot 
go back, i.e. the atom is caught between the diode and 
the trap wall. Note that the implementations of an atom 
diode presented in 0, [H-HI are designed in such a way that 
the internal state of the atom is changed during the pass- 
ing of the diode but it will be restored at the end in such 
a way that the atom is in the same internal state before 
and after crossing the diode. The diode, which is in prin- 
ciple a semi-penetrable barrier, behaves like a wall for the 
captured atom and continues moving to the left without 
changing its velocity. During this the atom is bouncing 
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off the slowly-moving diode and in such a way the ki- 
netic energy of the atom is reduced. Note that the total 
energy is not conserved during this process (a Hamilto- 
nian describing the system would be time-dependent). 
So the diode transports the caught atom to the bottom 
of the potential without increasing the kinetic energy of 
the atom, on the contrary, the kinetic energy of the atom 
might be even more decreased. This can be heuristically 
compared with a ball bouncing on a horizontal racket 
where the racket is slowly moved down in the gravity 
field: Then the ball is not bouncing "more" if the racket 
is arrived at the bottom. Note that even if the atom 
has reached the bottom of the trap, diode still continues 
moving to the left. Finally, the atom is in the same ex- 
ternal potential (this is different for example from [ll|) 
and in the same internal state as it was initially, but its 
total energy has been reduced. 

If the process acts on the whole ensemble of non- 
interacting atoms -of a given temperature initially- then 
the energy of every atom is reduced and therefore the 
ensemble is cooled. 

It is important to underline that -because of the irre- 
versible step- this cooling method is fundamentally dif- 
ferent from velocity reduction by collision with a moving 
wall which is not "real" cooling 12;]. 

Several variations of this cooling method have been 
proposed and applied in various experiments, both for 
cooling atoms p^ - iT7| and for cooling molecules [H, . 

As a variant to the scheme shown in Fig. [TJ a small 
optical dipole trap has been used in experiment flij . 
The small trap is constructed in the form of a square 
box [l4| and therefore we will simply call it "box" in 
this paper (not to be confused with the external poten- 
tial trap). The key idea is that this small box is moved 
slowly though the external potential. In an ideal setting, 
the walls of the box consist of atom diodes: if the atom 
has crossed a wall of the box and the kinetic energy of the 
atom is lower than the threshold energy of the "diodic" 
wall (called the threshold energy of the box in the fol- 
lowing) then the atom is irreversibly caught by the box. 
Note that in the ideal case the caught atom is in the same 
internal state as before the catching, see also In 
such a way, a similar effect as shown in Fig. [T] may be 
achieved. 

It is important that the box is moved in such a way 
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this trap is 



FIG. 1. Schematic representation of the process for a single 
atom: snapshots for different, increasing times (from left to 
right); the passing direction of the diode is from left to right. 




FIG. 2. (Color online) Schematic representation of the two 
different two dimensional traps considered: (a) wedge trap, 
(b) harmonic trap; the red-filled squares indicate the "box". 



that it traps a maximal number of atoms. While this 
optimal trajectory of the box is straightforward in a one- 
dimensional setting, this is not obvious in higher dimen- 
sions. In [20], a two-dimensional wedge trap has been 
examined. In that paper the box was assumed to be 
at rest and the optimal position of the box was deter- 
mined. The main goal of our paper is to extend this work 
[2fj| and show that by a slowly moving box the cooling 
efficiency can be improved compared to a box at rest. 
Therefore, we will examine a simplified model of catch- 
ing process in a two-dimensional wedge trap as well as in 
a two-dimensional harmonic trap, see Fig. [2J 

The paper is organized as follows: In Sect. [Til the 
model will be explained in detail. We will then examine 
the setting of a two-dimensional wedge trap in Sect. IIIII 
and of a two-dimensional harmonic trap in Sect. IIVI 



II. THE MODEL 

We examine two different two-dimensional trap geome- 
tries. The first one is an idealized two-dimensional, sym- 
metric wedge trap consisting of perfectly reflecting walls, 
see Fig. [5^. Gravity is acting in — y direction. The Hamil- 
tonian of a single atom in this trap is 



w 



1 
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(pl +pI) + mgy + V(x,y), 



(1) 



where V(x, y) = if y > 0, |a;| < y tan a and V(x, y) = oo 
otherwise. Moreover, m is the mass of the atom and g is 
the gravitational acceleration at the earth's surface. 

The second trap is a two-dimensional harmonic trap, 
see Fig. [2b- The Hamilton function of a single atom in 
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where uj is the frequency of the harmonic trap. 

In both cases, we approximate the motion of the atoms 
as a classical motion and in addition, we assume that the 
atoms inside the trap do not interact with each other. 
So the trajectory of a single atom can be calculated by 
solving Hamilton's equations 

dx dH dy dH dp x dH dp y dH 
dt dp x ' dt dp y dt dx' dt dy ' 

(3) 

The initial state of a single atom should be distributed 
concerning the canonical probability distribution with 
temperature T; 



Pi(x,y,p x ,p y ) = -^exp [ 



-/3H(x,y,p x ,p y )] 



(4) 



where Z is the canonical partition function, /3 — l/^^T 1 ;) 
and ks is the Boltzmann constant. 

Our simplified model of the catching process is similar 
to and motivated by [20|. We model the "box", as a 
square box in space with a width 2wb in x and in y 
direction, i.e. wb is the half width of the box. The box 
can move freely throughout the trap. We assume that the 
"diodic" walls of box are infinitely thin. Moreover, the 
"diodic" walls of the box have a certain threshold energy 
Eb such that their "diodic" behavior breaks down if the 
kinetic energy of the atom is above this threshold energy. 
If the atom has crossed -at some time- the walls of the 
box and the kinetic energy of the atom in the box rest 
frame is smaller than the threshold energy Eb then the 
atom is caught irreversibly by the box for all times. If 
the kinetic energy of the atom is too large (larger than 
Eb) then the "diodic" property of walls breaks down and 
the atom can escape the box such that it might have a 
second chance later to be caught. 

The algorithm is now the following: For a single nu- 
merical run, we choose the initial state of a single atom 
randomly with respect to the probability distribution Eq. 
([4]). Then we calculate numerically the evolution of the 
atom until a final time tt while always checking if the 
atom is in the region of the box during its motion. If 
this is the case, we are checking if the kinetic energy of 
the atom in the box rest frame is smaller than the box 
threshold energy Eb- If this is true, then the atom is 
caught and the run is finished, otherwise the motion of 
the atom continues and it may get trapped later. We 
repeat such a single run Ni = 10 6 times, get the to- 
tal number of trapped atoms Nb and finally the relative 
number of trapped atoms, i.e. the fraction F = Nb/Ni. 

The goal of this paper is to optimize this fraction F 
by varying the trajectory of the box in the trap. As we 
will see now, from this follows also an optimization of 
the cooling efficiency if the initial conditions and the box 
parameters are fixed. 
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One can define a cooling efficiency as the compression 
of the phase-space density. Phase-space density in the 
context of cooling can be defined as the number of atoms 
in a box with sides of one "thermal de Broglie wave- 
length" [21| . The "thermal de Broglie wavelength" A can 
be defined as 



A = 



(5) 



where m is the mass of the atom. The phase-space 
density in a two-dimensional setting is then defined as 
H = rtA 2 , where n is the spatial density of atoms. The 
relative change in phase-space density is therefore 



Hi Hi T f '' 



(6) 



where \ii is the initial and fj,f is the final phase-space 
density. Tj and Tt are the initial and final temperature, 
respectively. 

The initial distribution will be a canonical one with 
temperature Tj. The initial spatial density m is strictly 
speaking not constant in space. Nevertheless, we set ap- 
proximately m — Ni/Ai, where Ni is the initial number 
of atoms. To get Ai, we fix a small e > and define 
the initial region Ai such that the atom is initially in Ai 
with probability 1 — e. The final spatial density is set as 
nf = Nb/Ab, where Nb is the number of trapped atoms 
and Ab = 4u; 2 j is the area of the box. 

The cooling efficiency 77 is then defined as 



V = logi 



log 



Nb_ A^T^ 
Ni A B T f 



(7) 



Here Ei — ksTi. For simplicity, the final temperature 
is approximated by the threshold energy of the box, i.e. 
Tj = E B /k B - The cooling efficiency depends on the 
fraction of trapped atoms F = Ns/Ni, the ratio of the 
area of the thermal system to the area of the box, Ai /Ab , 
and the ratio of the initial to final energy, Ei/Es- 



V = log! 



,Ai_Ek_ 
A B E b 



(8) 



If the initial conditions and the box parameters are fixed, 
it will be therefore sufficient to optimize the fraction of 
trapped atoms F. The goal of the rest of the paper will 
be to optimize this fraction by choosing different box 
trajectories. 



III. BOX TRAJECTORIES IN A WEDGE TRAP 

First we study the wedge trap, see Fig. According 
to Eq. (j4]), the initial canonical distribution is given by 



Pi,w(x,y,Px,P y ) 



mg X[- 



y tan a, y tan 



4-7T tan a 

x exp {-(3 [( Px 2 + p y 2 )/(2m) + mgy] } , (9) 



where xj( x ) is the indicator function with Xj( x ) = 1 if 
x G J and Xj( x ) = if x (£ J. It is convenient to define 
a characteristic length I, a characteristic velocity v and 
a characteristic time r by 



I = 



k B Ti 
mg 



knT, 



l k B T t 
mg 2 



(10) 



In the rest of the paper we are assuming 87 Rb atoms, 
i.e. m = mass( 87 Rb) = 1.44316 x 10~ 25 kg. The initial 
temperature should be Tj = 100 /iK and g = 9.78 m/s 2 
(equator). The characteristic values are then 



I = 978 fim, 



9.78cm/s, r = 10ms. 



fill 



Box at rest 



The center of the box is moving with a trajectory 
ys(t)). As a reference case we first assume that 
the box center is placed at xb — at rest, i.e. its ve- 
locity is vb.x — v B,y = 0. Fig. [3^, shows the trapping 
fraction F versus different box positions ys for different 
combinations of wedge angle a and box half width Wb- 

The optimal height y op which gives the maximal frac- 
tion of trapped atoms for different box half widths wb is 
shown in Fig.[3jD. The error bars are defined by the range 
in which the maximal trapping fraction F decays by an 
amount of 1/ y/Ni, where Ni is the number of particles 
used in the numerical simulation as defined above. 



B. Box with linear motion 

We want to examine if a moving box can produce a 
higher fraction F of trapped atoms than the box at rest. 
First, we consider a linear motion of the box given by 



X B (t) = V B ,x{t 

y B (t) = v B , y (t 



*//2), 
tf/ 2 ) + Vop- 



(12) 



Here y op is chosen to be equal to the optimal value for the 
box at rest shown in Fig.[3]D (for the corresponding wedge 
angle and box half width). Because of mirror symmetry 
(x o —x) of the setting we can restrict to the case vb, x > 
0. 

Fig. [4] shows the resulting fraction F versus different 
box velocities for different wedge angles. The box half 
width is fixed at wb/1 — 0.35. The square symbols mark 
the maximal fraction F which can be achieved with a box 
at rest while the dots mark the maximal fraction which 
can be achieved with a box moving with the trajectory 
Eq. (jT^j) . It is clearly seen that the linear moving box 
can capture a larger fraction F of atoms than the box at 
rest (for all the three different wedge angles a = 30°, 45° 
and 60°). 

The optimal velocity in the first case a = 30° (Fig. |4ji) 

is 



vb,x/v 

V B ,y/v 



0.06 
0.12 



0.13 sin 30° . 
0.13 cos 30° 



(13) 
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FIG. 3. (Color online) Wedge trap with the box at rest at 
xb = 0. (a) Fraction F of trapped atoms versus box positions 
yB'- a — 30°, wb/1 = 0.1 (solid blue line); a = 45°,wb/1 = 
0.1 (dashed red line) ; a = 60°, wb/1 = 0.1 (dotted green line); 
a = 30°, wb/1 = 0.35 (thick solid blue line); a = 45°, wb/1 = 

0. 35 (thick dashed red line); a = 60°, wb/1 = 0.35 (thick 
dotted green line), (b) Optimal position y op of the box at 
rest for different wedge angles (blue, red resp. green plus 
signs connected with lines), the black dots correspond to the 
straight line y op = u>b(1 + tan a)/ tan a, see text for more 
details. tf/r = 20, E B /Ei = 0.1. 

In the case a = 45° (Fig. 0b), we get for the optimal 
velocity 

vb.x/v = 0.08 PS 0.11 sin 45°, . . 

v B , y /v = 0.08 « 0.11cos45°. [ > 

It is important to notice from Eqs. (p~3|) and (IT4|) that 
we get in both cases the maximal fraction if the box is 
approximately moving parallel to one wedge side. 

This is different in the case a = 60° (Fig. gfc). The 
optimal velocities are now 

v B ,x/v = 0.04 « 0.072 sin34°, , , 

v B \ y /v = 0.06 w 0.072 cos 34°, ^ b > 

1. e. the box motion is not parallel to the wedge side in 
this case. Nevertheless, the linear moving box traps in 
all cases more atoms than the corresponding box at rest. 

Therefore, the box moving with the linear trajectory 
enhances the trapping fraction than the box at rest. In 
the following, we will study the linear moving box in more 
detail to see if this result remains true for other box half 
widths. 

We first look at the cases of a wedge angle a — 30° and 
a = 45°, respectively. Figs. [6] and [7] show the fraction F 
of trapped atoms versus different box half widths wb- 
Figs. [5^. and [7^, correspond to the final time tf/r = 20 
while Figs.JBJj and[7)3 correspond to tf/r — 40. 




(c) 

FIG. 4. Wedge trap. Box moving with the linear trajectory 
Eq. (|12[1 . Fraction F of trapped atoms versus box velocities: 
(a) wedge angle a = 30°, y op /l = 0.7; (b) wedge angle a = 
45°, yop/l = 0.6; (c) wedge angle a = 60°, y op /l = 0.75. 
tf/r = 20, E B /Ei = 0.1, wb/1 = 0.35. The square symbols 
mark the maximal fraction for a box at rest, the dots mark 
the maximal fraction for a linear moving box. 




FIG. 5. (Color online) Schematic representation of the box 
trajectories Eqs. (|16[) and (JTTJ) , the initial position of the box 
is shown, see text for further details. 

The results for a box at rest {xb =0, vb,x = ^B,y = 0) 
are shown as a reference case for tf/r= 20 as well as for 
tf /t = 40 (plus signs connected by thick blue line) with 
the box coordinate yB — Hop and y op shown in Fig. (3Jd 
(which is optimal for tf/r = 20). It can be seen in Figs.[S] 
and [7] that for a box at rest the fraction does not change 
significantly if the final time tj increases. 

Motivated by the previous results, we are again con- 
sider a box which moves linearly in direction of the right 
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wedge side, i.e. its velocity is vb, x — vsma,VB,y — 
v cos a. The box trajectory should cross the y-axis at 
y op , where y op is again the position shown in Fig. [3]o. 
Nevertheless, the box should now start at time t = 
directly outside the wedge trap, see Fig. [5] Note that 
therefore the y axis is no longer crossed at time tf/2 by 
the box center. The resulting trajectory is 

Xs(t) = v sin at - ^ (w B + (w B + y op ) tana) , , „, 
UB{t) = v cos a t — 1 (wb — Vop + wb cot a) . 

The results for such a moving box for tf/r = 20 as well 
as for tf/r — 40 are shown in Figs. |B][7] (crosses connected 
by red dashed line) . It can be seen that the fraction F is 
-in all cases- much larger for such a linearly moving box 
compared to the box at rest. 

However, these box trajectories still depend on the val- 
ues y op which have to be obtained by numerical optimiza- 
tion. The goal is now to find an analytical approximation 
for y op . Therefore, we are considering a box moving par- 
allel to the right wedge side while touching the wedge 
side with its right, lower corner (see also Fig. [5]). Then 
we get the analytical value y op = wb(1 + tana)/ tana. 
Note that these values for y op are also plotted in Fig. [3b 
and we can see that this is also a rough approximation 
for the numerical determined optimal position y op . This 
"analytical" box trajectory is now 



XB{t) = w sin at — wg(l 
= vcosat. 



- tana) 



(17) 



The resulting fraction F using the box trajectory Eq. [T7] 
can be seen in Figs. |BE] (dots connected by thick green 
dotted line) . We find that the box traps even more atoms 
than with the trajectory Eq. (fl6|) considered in the previ- 
ous paragraph and it has the advantage that no numerical 
determined value of y op is required. 

For completeness, the fraction F versus the box half 
widths for a = 60° is shown in Fig. [8] The results are 
shown for a box at rest (plus signs connected by a thick 
blue line) and a box moving linearly with trajectory Eq. 
(QU as well as Eq. ([TT]). The value of y op in Eq. ([T6]) is 
chosen from Fig.[3j3. We see that the moving box catches 
more atoms than the box at rest. In contrast to the cases 
a = 30° and a = 45°, the analytical box trajectory Eq. 
(flT]) is here a good choice only for tf/r = 40. 



C. Box moving with wriggled trajectory 

We want to underline that a linear trajectory might 
not be optimal. However, the advantage of a linear mo- 
tion of the box is that it might be easier to implement 
experimentally a linear motion than a more complicated 
motion of the box. 

Nevertheless, as an example, we shall also try a differ- 
ent trajectory of the box, a wriggled one, given by 



%B{t) = y b (t) tan a cos(u>wt) , 
ys(t) = ywo + [wb - ywo] jj ■ 
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FIG. 6. (Color online) Wedge trap, a = 30°. Fraction F of 
trapped atoms versus box half widths wb'- (a) tf/r = 20, (b) 
tf/r = 40; different box trajectories: box at rest (plus signs 
connected by thick blue line), box moving with trajectory Eq. 
(|16[) (crosses connected by red dashed line), box moving with 
analytical trajectory Eq. (|17[) (dots connected by thick green 
dotted line), box moving with wriggle trajectory Eq. (|18[) 
(circles connected by black solid line), v/u — 0.13, ywo /I = 
2.0, uj w t = 0.25, E B /Ei = 0.1. 




0.2 0.25 0.3 0.35 
Wr/I 



(a) 




0.35 



(18) 



(b) 



FIG. 7. (Color online) Wedge trap, a = 45°. Fraction F of 
trapped atoms versus box half widths wb'- (a) tf/r = 20, (b) 
tf/r = 40. v/v = 0.11, ywo/l = 1.5, ujwt = 0.2; see Fig. E] 
for more details. 
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A. Box at rest 

As a reference, we first look again at the case that the 
box center is placed at xb = at rest, i.e. Vb, x — v B,y — 
0. The box threshold energy is Eb/Ei — 0.1 and the 
final time is tf/r — 60. The box coordinate ys = y op 
is numerically optimized such that the resulting fraction 
of trapped atoms F is maximal. Fig. [9^ shows the trap- 
ping fraction F for different box half widths wb (plus 
signs connected by thick blue line). The optimized box 
coordinate y op is shown in Fig. [9}) (plus signs connected 
by thick blue line). The error bars are defined by the 
range in which the maximal trapping fraction F decays 
by an amount of 1/ ^/Ni, where JVj is the number of par- 
ticles used in the numerical simulation. From the error 
bars, we can see that for small wb the result is relatively 
independent of the exact box position ys- 



B. Box with linear motion 



FIG. 8. (Color online) Wedge trap, a = 60°. Fraction F of 
trapped atoms versus box half widths wb- (a) tf /t = 20, (b) 
tf/ T = 40. v/v = 0.11, ywo/l = 1-5, lj w t = 0.2; see Fig. E] 
for more details. 

The box center starts at the right wall of the trap, i.e. 
x B (0) — ywo tan a, ?/b(0) = ywo, is then moving down 
wriggling and ends at y(tp) — wb, i.e. the box half 
width. The resulting trapping fractions F can also be 
seen in Figs. EMU (circles connected by black solid line). 
The parameters for the box trajectory have been chosen 
such that the velocity of the box is limited v/v < 0.265 
at all times. In the cases a = 30° and a = 45° we get an 
increased fraction of trapped atoms for i//r = 40. 

IV. OPTIMIZING BOX TRAJECTORIES FOR A 
HARMONIC TRAP 

Now, we study the harmonic trap with frequencies 
uj x = uj y = Li, see Fig. ^jp. The initial state of an atom is 
chosen concerning the canonical distribution 

( \ 
Pi,H(x,y : p x ,p y ) = 

x exp {-p [{pi +pg)/(2m) + mu: 2 (x 2 + y 2 )/2] } (.19) 

Again, it is convenient to define a characteristic length /, 
a characteristic velocity v and a characteristic time r by 

l= l -J^, » = J^, r = i. (20) 
cu V m V m to 

In the rest of the paper we are again assuming Rb 
atoms, the initial temperature shall be Tj = 100 fiK and 
uj = 50 x 2ir Hz. The characteristic parameters are then 

Z = 311/im, ^ = 9.78cm/s, r = 3.18 ms. (21) 



We want to examine if a linearly moving box can pro- 
duce a higher fraction F of trapped atoms in the case 
of a harmonic trap. Therefore, we are now looking at a 
box moving in such a way that its trajectory crosses the 
point x = 0, y = y c at time tf/2. Because of symmetry, 
we can restrict to the case vb,x > 0, vb, v = (if we ne- 
glect the rotation of the box itself). We consider a linear 
trajectory of the moving box as 

x B {t) = v B , x (t-tf/2), , 22 > 

VB{t) = Ve- 
in Fig. [TUJ the resulting fraction F is shown for the box 
half width wb/1 = 0.2. The box symbol marks the max- 
imal fraction for a box at rest while the dot marks the 
maximal fraction for a linearly moving box. It can be 
seen that the linearly moving box traps more atoms than 
the box at rest. 

We also calculated the fraction for different box half 
widths, the result can be seen in Fig. [S^. For every 
box half width, the position y c and the velocity v B ,x 
have been optimized and the optimal values are shown 
in Fig. [SJd and Fig. [SJ: (crosses connected by red dashed 
line). The error bars are defined in the same way as ex- 
plained in the last subsection. In Fig. [9^ it can be seen 
that in all cases the moving box traps significantly more 
atoms than the box at rest. Good approximations for the 
optimal position y c and the velocity vb,x are 

y c /l - 0.55, (23) 
v B>x /v = 0.025 + Q.2hw B /l, (24) 

which are also shown in Fig.[U> and Fig. (black dashed 
line). These approximations lead to an analytical box 
trajectory of 

x B (t) = (0.025 v + 0.25w B /r)(i-t//2), . . 
y B (t) = 0.55 1 . [Zb) 
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FIG. 9. (Color online) Harmonic trap, data versus different 
box half widths wb- (a) Fraction F of trapped atoms, box at 
rest with optimized y op (plus signs connected by thick blue 
line), linear moving box with optimized vb,x and optimized y c 
(crosses connected by red solid line), linear moving box with 
trajectory Eq. (|25|) (dots connected by thick green dotted line, 
on top of red line), box moving with a helix trajectory (circles 
connected by black solid line), (b) Optimized box parameter 
y op resp. y c (crossing point of y-axis): for a box at rest (plus 
signs connected by thick blue line), for a linearly moving box 
(crosses connected by red dashed line), approximation Eq. 
(|23[l (black dashed line), (c) Optimized box velocity vb,x'- 
for a linearly moving box (crosses connected by red dashed 
line), approximation Eq. (I24[) (black dashed line). Es/Ei = 
0.1, t f /r = 60. 



The resulting fraction F for a box moving with trajec- 
tory Eq. (j2"5l) is also shown in Fig. [3k. (dots connected by 
thick green dotted line on top of the red line), it is in- 
distinguishable from the result obtained with trajectory 
Eq. (I22[) with optimized parameters y c and «b. re- 
in the following we will show that this trajectory Eq. 
(1231) is also a very good choice for different box threshold 
energies Eb- The fraction F versus different box thresh- 
old energies Eb for the box at rest (plus signs connected 
by thick blue line) and the box moving with trajectory 
Eq. (f2"5j) (dots connected by thick green dotted line) are 
shown for the box half width wb/1 = 0.2 in Fig. [TTa and 
for wb/1 — 0.35 in Fig. ITTb . It can be seen that even for 
other box threshold energies Eb the linearly moving box 
traps more atoms in the same time than the box at rest. 



0.1 - 




FIG. 10. Harmonic trap, box moving with the linear tra- 
jectory Eq. (122[) , fraction F for different positions y c and 
box velocities vb,x', vb, v = 0; wb/1 = 0.2; the square 
symbol marks the maximal fraction for a box at rest, the 
dot marks the maximal fraction for a linear moving box. 
E B /E l = 0.1, t f /T = 60. 
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(a) 




EJE, 



(b) 



FIG. 11. (Color online) Harmonic trap, fraction F versus 
box threshold energies Eb: box at rest (plus signs connected 
by thick blue line), box moving with the linear trajectory 
Eq. (|25|l (dots connected by thick green dotted line), box 
moving with a helix trajectory (circles connected by black 
solid line), (a) Box half width Wb/1 = 0.2, (b) box half width 
wb/1 = 0.35. tf/r = 60. 



C. Helix trajectories 

Again, we want to emphasize that the linear trajec- 
tory above might not be optimal but it is probably easier 
to implement experimentally than a more complicated 
motion. 

Just as an example, we try a different, more compli- 
cated box trajectory. The box is moving with a helix 
trajectory given by 



= x H (l~t/t f )cos(uj H t), 
Vsit) = x H (l-t/t f )shx(u! H t), 



(26) 
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where xh and ujh are the key parameters of the box 
trajectory. The results can be seen in Fig. [!|] as well as 
in Fig. [TT] (circles connected by black solid line). The 
parameters of the helix trajectory are xh/1 = 1-9 and 
ujh = 0.1. For larger box half width wb (see Fig. [Tib) 
this trajectory helps to get a larger fraction of trapped 
atoms while there is no improvement for smaller box half 
width. 



V. SUMMARY AND OUTLOOK 

We have examined a model for cooling in a two- 
dimensional wedge trap and a two-dimensional harmonic 
trap. During the cooling procedure, the atoms are cap- 
tured by a small area surrounded by "diodic" walls, called 
"box" , which moves through the wedge trap and the har- 
monic trap, respectively. We have examined different 
box trajectories with the goal to maximize the fraction 



of trapped atoms, this leads also to an increased cooling 
efficiency. We have shown that the fraction of caught 
atoms can be increased using a moving box compared 
to a box at rest which was examined in an earlier work 
[20( | . We have also optimized the parameters of the box 
trajectory where we restricted ourselves mainly to lin- 
ear box motions due to its possible easier experimental 
implementation. In a future work, we will consider the 
question of the optimal general box trajectory in more 
detail and also taking quantum effects into account. 
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